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Abstract 

The evaluation of a tunneling tail by the Herman-Kluk method, which is a quasiclassical way to 
compute quantum dynamics, is examined by asymptotic analysis. In the shallower part of the tail, 
as well as in the classically allowed region, it is shown that the leading terms of semiclassical eval- 
uations of quantum theory and the Herman-Kluk formula agree, which is known as an asymptotic 
equivalence. In the deeper part, it is shown that the asymptotic equivalence breaks down, due to 
the emergence of unusual "tunneling trajectory", which is an artifact of the Herman-Kluk method. 

PACS numbers: 03.65.Sq, 05.45.Mt, 02.70.-c 



*Electronic address: tanaka@phys.metro-u.ac.jp 



1 



Even nowadays, it is still impossible to carry out serious numerical investigation of quan- 
tum dynamics, even for modest (e.g., 10-100) degrees of freedom systems, with the present 
state of the art of computational technology, unless we take drastic approximations that 
need to be based on good physical insights. As a starting point to invent such a method, it 
is often employed the semiclassical approximation (i.e. asymptotic evaluation) of the path 
integral representation of time evolution operator (Feynman kernel) |l|. The semiclassical 
approximation, however, has difficulties due to the exponential proliferation of contributing 
trajectories and caustics (see, e.g. Ref. Jsj]) and due to Stokes phenomenon that requires to 
remove non-contributing complex- valued trajectories 0, even in few degrees of freedom sys- 
tems, when the corresponding classical system is chaotic. Furthermore, in the semiclassical 
method, the boundary conditions of the contributing classical trajectories both for initial 
and finial times, is troublesome (known as a root-search problem 0,O]), in computations for 
realistic systems such as atoms and molecules. In order to avoid the root-search problem, 
initial value representations (IVRs) of Feynman kernels are proposed. The IVRs imposes 
only the initial conditions on the classical trajectories. In Ref. Q, the earliest version of 
IVR is introduced by a change of integral variables to semiclassical Feynman kernel. A 
general framework of IVR is proposed by Kay p , who discussed various integral expressions 
(IEs) of approximate Feynman kernel, which are composed by classical trajectories that are 
emitted from real- valued initial conditions. The important guiding principle of Kay's theory 
is that the leading semiclassical expressions of both an IE and exact quantum theory must 
agree. This is called the asymptotic equivalence jflHl. Furthermore, Kay argued that several 
known IVRs (IEs), including thawed Gaussian approximation [sj], cellular dynamics [9f, and 
the Herman-Kluk (HK) formula ^(|, are asymptotically equivalent with quantum theory. 
Nowadays, a lot of numerical investigations of quantum dynamics, including rather realistic 
systems, employ IVRs, in particular, the Herman-Kluk method |l l| . 

The limitation of IVRs, however, is not clear in particular, in the descriptions of 
classically forbidden processes, e.g., tunneling processes, whose conventional semiclassical 
treatments need to take into account the contributions from complex-valued classical tra- 
jectories. Numerical experiments to reproduce tunneling tails by IVR approaches suggest 
that the "shallow" side of tunneling tails is tractable [13]. On the other hand, concerning to 
the "deep" side, the IVR approaches have difficulties. Kay's semiclassical analysis of 0{h 2 ) 
error term of Herman-Kluk method reveals that the magnitude of the error is controlled 
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by the complex-valued classical trajectories However, these works are not conclusive. 
First, since the tunneling tails are exponentially small, the corresponding error analysis also 
requires to treat exponentially small errors. Hence, the 0(h 2 ) error term, which is satisfac- 
tory in classically allowed region, is too large. Second, the border between the shallow side 
and the deep side of the tunneling tail has been unknown. In order to clarify the limitation 
of IVR, the identification of the border is inevitable. We remind that asymptotic (i.e., semi- 
classical) analysis has an ability to treat exponentially small quantities. This suggests some 
asymptotic approach may reveal the limitation of IVRs with much better accuracy. 

In this paper, we examine an evaluation of a tunneling tail, by the Herman-Kluk formula, 
with asymptotic analysis. Here we focus on the Feynman kernel, rather than energy spectra 
or correlation functions. This facilitates to identify the origin of discrepancies. 

We here examine a single degree of freedom system that is described by a Hamiltonian 
H = —gp 3 /3, where q and p are the position and the momentum of the system, respectively. 
We assume that the strength of folding g is positive. This is a canonical model that describes 
a nonlinear folding process in classical phase space (see, Fig.^(a)) 3]. When both the initial 
and the final states of Feynman kernel are eigenstates of the position operator, the nonlin- 
ear folding dynamics produces a caustic. Note that in generic, nonlinear systems induce 
foldings in general. Dynamics locally around each foldings are described by the canonical 
Hamiltonian with appropriate rotations and scaling in phase space (see, e.g., Fig. d(b)). 

The Feynman kernel of a time evolution of interval [0, r] (r > 0) in the position repre- 
sentation K(q) = (q\ exp(—iHr/h)\q = 0) is expressed exactly with Airy function 



where I = (ft^gr) 1 / 3 is a characteristic length for a penetration into the classically forbidden 
region q > 0. The asymptotic evaluation of the integral representation of the Feynman 
kernel 



gives a leading semiclassical approximation K^c, which are composed by classical trajec- 
tories. On one hand, in the classically allowed region q < 0, we have a superposition of 
incoming and outgoing waves: 



K(q) = Mq/l)/l 
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FIG. 1: Evolutions of classical manifolds in phase space, (a) With the folding Hamiltonian 
H = —gp 3 /3 j3(| . The initial manifold (thin line) is at q = 0. After a time interval r, the manifold 
folds (thick line). The caustic in the position representation is at q = 0. (b) With a nonlinear 
oscillator H = \p 2 + V(q), where V(q) = D{{1 - e'^) 2 - 1} + |(1 - e)q 2 , e = 0.975, A = l/Vu, 
and D = e/(2A 2 ) (Contour lines of V(q) are thin) JJ|. The initial manifold {(q,p)\q = 9,-3 < 
p < 3} (vertical line) mimics an eigenstate of the position operator, with an energy cutoff. The 
corresponding final manifold (thick line), at t = 18, has two prominent caustics in the position 
representation . 

The corresponding classical trajectories are characterized by their momentum, which are 



conserved quantities: p±(q) = ±^/\q\/(rg). On the other hand, in the classically forbidden 
region q > 0, we have a tunneling tail: 

1 



Ksc(q) 



2 v /^(g//)V4 



exp 
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The momentum of the tunneling trajectory is pure imaginary: po(q) = 1 qj (rg). Between 



these regions, at q = 0, these classical trajectories merge to produce a caustic. The change 
of asymptotic expansions for different signs of q is controlled by Stokes phenomenon 
In the following, the corresponding Herman-Kluk kernel 10( is examined: 
'HKv \ ... f f d( lodpo 



K = J J ^t^^^,Pr))C{q Q ^T) 

xe iS ^ p ^ h (^(q , Po )\q = 0), (5) 

where 7 (> 0) determines the width of Gaussian packet (q\^(qo,Po)) = 
(27/7r) 1 / 4 exp{— 7(g — q ) 2 + ipo(q — qo)/h}i (Qt,Pt) is the classical trajectory at time t — r, 
emitted from (qo,po) at t = 0, S T (qo,po) is the classical action along the time evolution, and 
C(qo,Po,T~) = {dq T /dq + dp T /dp — 2ih'ydq T /dpo — (2ih'y)~" 1 dp T /dqo} 1 / 2 /\/2. For the folding 
Hamiltonian, C(qo,po,r) have a branch point at p\ = i/{2fvyTg). After the integration of 
the variable go, K nK (q) (jSJ) becomes 

K HK (q) = J^C(P^- MP \ (6) 
where C(p,r) = (1 — p/pi) 1 ^ 2 and 4> T {p) = l(q + T 9P 2 ) 2 /^ ~ ipq/h — irgp 3 /(3h). The 



integral © has three saddle points p = ±y^—q/(rg) and p\. The former momenta p = 
±^—q/(rg) correspond to the classical momenta p±(q) in the classically allowed region 
and po{q) in the classically forbidden region. The latter momentum pi, which is the branch 
point of C(p,r), appears only in the semiclassical analysis of K nK (q). Note that all the 
saddle points need not to make contributions to the semiclassical kernel, due to the Stokes 
phenomena. Actually, in the classically allowed region q > 0, there are the contributions 
only fromp = p±(q) (FIG.|2fa)). Hence the leading semiclassical evaluation of K (q) agree 
with Ksc{q)- Thus the asymptotic equivalence between Herman-Kluk kernel and quantum 
theory holds for q < 0]. 

In the semiclassical evaluation of K KK (q) in the classically forbidden region q > 0, there 
is a length scale Z 7 = 7~ 2 Z~ 3 /4 that divides the tunneling tail into two regions: a "shallow" 
region < q < Z 7 and a "deep" region q > Z 7 . In the shallow region, the semiclassical 
evaluation of K nK (q) has only a single contribution from the classical trajectory p = po(q) 
(FIG.Efb)). Hence the asymptotic equivalence between Herman-Kluk kernel and quantum 
theory holds both classically allowed region and the shallow tunneling region q < / 7 . This is 
a promising evidence that collections of real classical trajectories can describe shallow tun- 
neling tails, through IVR techniques. This is the first example to show that IVR technique 
can describe a classically forbidden process with an analytical argument. 
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FIG. 2: Locations of saddle points (closed and open circles) and integration paths (thick lines) 
for the asymptotic evaluations of the integral (JBJ), for (a) the classical region q < 0, (b) the shallow 
tunneling region < q < l~, and, (c) the deep tunneling region q > L. Dashed lines emanating from 
the branch point pi are branch cuts. Closed and open circles are contributing and non-contributing 
saddle points, respectively. 

On the other hand, in the "deep" tunneling region, we found a discrepancy: Due to 
Stokes phenomenon, the contribution from "conventional" tunneling trajectory p = po(q) 
disappears, and in turn, the contribution from the classical trajectory p = pi, which is 
an artifact of the Herman-Kluk kernel, appears (FIG. I2£c)). The resultant semiclassical 
evaluation of K KK (q) is 



I 



3/4 



r(3/4) 

2tt/( 7 / 2 ) 1 /4 V<?-/ 7 . 



2 v. ■ -v 3 "7 f ' ( 7 ) 

Note that q = Z 7 is a caustic, which is an reminiscent of the Stokes phenomena. At the same 
time, the asymptotic form of tunneling tail ~ exp(— 7g 2 /2) for q ^> Z 7 , which sensitively 
depends on 7, is qualitatively different from Ksc(q) ~ exp{— 2(g//) 3//2 /3} (JU). Thus the 
breakdown of the asymptotic equivalence is evident. 

In the argument above, the discrepancy between Herman-Kluk kernel and Feynman kernel 
comes from (1) a nonlinear folding dynamics in the corresponding classical phase space, and 
(2) the appearance of artificial tunneling trajectory. When the folding dynamics is not 
significant (this is the case before Ehrenfest time), there is a workaround to remove the 
contribution from the artificial tunneling trajectories, by adjusting the value of 7 in the 
Herman-Kluk kernel to push p\ deeper in the complex plane (see., Fig. |2J). Indeed, this is a 
known strategy, which is proposed by Kay, to reduce the magnitude of the error of Herman- 



Kluk method, and the strategy succeeds to a certain extent [1J, |17|. However, in generic, 
nonlinear systems, it will become difficult to carry out such workarounds in practice, due to 
the emergence of multiple foldings (see, e.g., Fig. [U(b)). 

We summarize this paper. With an exactly solvable model that describes nonlinear 
folding process in corresponding classical phase-space dynamics, we identified a boundary 
between a shallow and deep tunneling regions for the Herman-Kluk kernel. In the former 
region, the Herman-Kluk kernel and quantum theory are asymptotically equivalent. Hence 
there remains a hope that Herman-Kluk kernel describes classically forbidden process to a 
certain extent. However, in the deep region, the breakdown of the asymptotic equivalent is 
shown. Besides Herman-Kluk method, other IVR approaches, in particular, which are based 
on Kay's framework p, will have similar scenario on successes and failures in descriptions 
of classically forbidden phenomena. We expect that the nonlinear folding Hamiltonian em- 
ployed here is a good test-case not only for the Herman-Kluk method, but also for various 
semiclassical method. We remind that this model was employed to elucidate the limitation 
of single trajectory approximation of semiclassical coherent-state path integrals . At 

the same time, it is highly desirable to develop a convenient way to find the boundary be- 
tween the shallow and the deep regions of Herman-Kluk method, though the present analysis 
requires a full knowledge of the Stokes phenomena to determine the boundary. 
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FIG. 3: Comparisons of exact, Herman-Kluk (HK), and semiclassical Herman-Kluk (HKSC) 
evaluations of the Feynman kernel K(q), in particular, its tunneling tail q > 0. The penetration 
depth is I = 1. (a) With Z 7 = 1, all three theories are shown. At the (conventional) turning 
point q = 0, the semiclassical Herman-Kluk encounters caustic. In the shallow region < q < Z 7 , 
the discrepancy between quantum theory and Herman-Kluk method is not significant. On the 
other hand, at q = L, the semiclassical Herman-Kluk encounters another caustic. Note that the 
conventional semiclassical theory meets the caustic at q = 0, the conventional turning point, only. 
In the deep region q > Z 7 , the tunneling tails of quantum theory and Herman-Kluk formula take 
qualitatively different shape (see, Eq. © and Eq. 0). (b) Comparison of quantum theory and 
Herman-Kluk formula, with several values of £ 7 . In order to reproduce the shape of deep tunneling 
tail by Herman-Kluk method, we need larger Z 7 l^J ]. 
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